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Abstract 


Two 20-day, continental midlatitude cases are simulated with a three-dimensional 
(3D) cloud-resolving model (CRM) and compared to Atmospheric Radiation 
Measurement (ARM) data. This evaluation of long-term cloud-resolving model 
simulations focuses on the evaluation of clouds and surface fluxes. 

All numerical experiments, as compared to observations, simulate surface 
precipitation well but over-predict clouds, especially in the upper troposphere. The 
sensitivity of cloud properties to dimensionality and other factors is studied to isolate the 
origins of the over prediction of clouds. Due to the difference in buoyancy damping 
between 2D and 3D models, surface precipitation fluctuates rapidly with time, and 
spurious dehumidification occurs near the tropopause in the 2D CRM. 

Surface fluxes from a land data assimilation system are compared with ARM 
observations. They are used in place of the ARM surface fluxes to test the sensitivity of 
simulated clouds to surface fluxes. Summertime simulations show that surface fluxes 
from the assimilation system bring about a better simulation of diurnal cloud variation in 
the lower troposphere. 
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1. Introduction 


The representation of clouds is one of the most important uncertainties in general 
circulation models (GCMs) in simulating the global water and energy cycle (e.g., Cess et 
al. 1990; Zhang et al. 2005). Currently, cloud-resolving models (CRMs) are being 
incorporated into large-scale dynamic frameworks (e.g., GCMs) to facilitate the 
interaction between clouds and large-scale circulations in place of conventional cumulus 
parameterization (Emanuel and Raymond 1993) as an alternative approach (e.g., 
Grabowski 2001; Khairoutdinov and Randall 2001; Raymond and Zeng 2005; 
Khairoutdinov et al. 2005; Chem et al. 2005). This approach is referred to as a 
superparameterization or multi-scale modeling framework (MMF, Randall et al. 2003). 
The approach has merit, since no assumption is introduced on the causality between 
cumulus clouds and large-scale circulations. A key question is whether current CRMs can 
function in a MMF as expected. This question can be addressed by evaluating long-term 
CRM simulations against observations. Here long-term simulation means a period 
comparable to the long timescale in the approach to radiative-convective equilibrium 
(Tompkins and Craig 1998), which is around three weeks. 

Consider a CRM in an ideal MMF with no computational limits. The CRM can be 
represented in a computationally expensive way [e.g., a three-dimensional (3D) 
framework with sufficient gridpoints]. When the CRM is driven with prescribed large- 
scale forcing, the difference between the modeling results and observations is attributed 
to model physics instead of computational issues [e.g., a two-dimensional (2D) 
framework, insufficient gridpoints]. Thus, the difference between the modeling results 
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and observations provides insights on improving the CRMs that are currently being used 
in MMFs. In this direction, the Goddard Cumulus Ensemble (GCE) model is evaluated 
with data collected in the Atmospheric Radiation Measurement (ARM) program. Two 
20-day, continental midlatitude cases were selected for this purpose. 

CRM evaluation can be traced back two decades. Although real clouds and cloud 
systems are 3D, most CRMs used today are still 2D due to computer resources (Krueger 
1988; Xu and Randall 1996; Wu et al. 1998; and many others). Only a few 3D CRMs 
(e.g., Tao and Soong 1986; Lipps and Hemler 1986) have been used to study the response 
of clouds to large-scale forcing. Previous studies showed that the collective 
thermodynamic feedback effect, and the vertical transports of mass, sensible heat, and 
moisture were quite similar between 2D and 3D simulations (e.g., Tao et al. 
1987). Recently, several 3D CRM experiments were performed for 7-day periods for 
tropical cloud systems with large horizontal domains (500x500 km 2 ) at the National 
Center for Atmospheric Research (Grabowski et al. 1998; Wu et al. 1998), NOAA/GFDL 
(Donner et al. 1999), Colorado State University (Khairoutdinov and Randall 2003) and 
NASA Goddard Space Flight Center (Tao 2003). Grabowski et al. (1998) found that 
cloud statistics as well as surface precipitation are significantly different between 2D and 
3D simulations when clouds and radiation are fully interactive. At Goddard, Tao (2003) 
found that the 3D-modeled water vapor budget is in much better agreement with 
observations in the lower troposphere than its 2D counterpart. 

In this study, the 3D GCE model is used to simulate continental midlatitude clouds for 
twenty days, longer than most previous 3D CRM simulations, to evaluate the simulated 
cloud residue (i.e., cloud ice, cloud water and water vapor). Since cloud residue is a 
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natural consequence of cloud growth, temporal and accumulated model surface 
precipitation is compared with observations, because surface precipitation rate can be 
regarded as a measure of cloud growth. Special attention is paid to the sensitivity of 
clouds to dimensionality. 

Model evaluation depends strongly on observational data quality (Moncrieff et al. 
1997). Some recent field experiments [e.g., the ARM Spring 2000 Intensive Operational 
Period (IOP) and the ARM 2002, International H 2 O Project (IHOP) in 2002] have 
provided comprehensive observational data (e.g., Ackerman and Stokes 2003; 
Weckwerth et al. 2004) for model evaluation. The observed cloud properties from the 
ARM Spring 2000 IOP, for example, are the best in the fourteen-year history of the ARM 
program (Xu et al. 2005). 

This study involves two twenty-day observation periods during two field experiments. 
One is from the ARM IOP 2000 and the other IHOP 2002. Some of the cloud systems in 
these experiments have already been studied for specific purposes, such as the life cycle 
of convective clouds and as a setting for a model inter-comparison (Wakimoto et al. 
2004; Xie et al. 2005; Xu et al. 2005). Xie et al. (2005) and Xu et al. (2005) focused on 
two short periods in the ARM IOP 2000 to compare four 2D CRMs as well as eight 
single-column models with observations. They suggested that 3D CRM simulations be 
done to narrow down the origins of differences between the 2D CRMs and observations. 

This paper evaluates clouds in long-term CRM simulations. It is organized as follows. 
Section 2 describes the models and observational data used. Sections 3 and 4 analyze the 
modeling results for springtime and summertime cloud systems, respectively. And, 
sections 5 and 6 give a discussion and summary, respectively. 
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2. Design of numerical experiments 

a. Model description 

A single-column CRM, which differs from the two-column models designed to 
account for the interaction between convection and large-scale circulations (Nilsson and 
Emanuel 1999; Sobel and Bretherton 2000; Raymond and Zeng 2000, 2005), is used here 
to test the response of clouds to prescribed large-scale forcing derived from observational 
data. The model structure follows previous ones (e.g., Johnson et al. 2002; Xie et al. 
2005; Xu et al. 2005) except for surface flux input. In the current framework, clouds are 
simulated with the GCE model, large-scale forcing data come from observations, and 
surface fluxes in the lower boundary come from either observations or a land data 
assimilation system. 

The GCE model is detailed in Tao and Simpson (1993) and Tao et al. (2003), which 
describes its development and main features. Its application to studies of precipitation 
processes and improving satellite retrievals can be found in Simpson and Tao (1993) and 
Tao (2003). The model is non-hydrostatic and anelastic. It can be used in two- or three- 
dimensions with cyclic lateral boundary conditions. Solar and infrared radiative transfer 
processes (four-stream discrete-ordinate scattering) are included. Their impact on cloud 
development associated with cloud-radiation interaction has been assessed (Tao et al. 
1996). Subgrid-scale (turbulent) processes in the model are parameterized using a scheme 
based on Klemp and Wilhelmson (1978) and Soong and Ogura (1980). The effects of 
both dry and moist processes on the generation of subgrid-scale kinetic energy have been 
incorporated. A three-class ice formulation (3ICE), namely that by Lin et al. (1983), was 
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used. The sedimentation of ice crystals was recently included in the GCE based on 
Heymsfield and Dormer (1990) and Heymsfield and Iaquinta (2000) and was discussed in 
detail in Hong et al. (2004). All scalar variables (temperature, water vapor, and all 
hydrometeors) are calculated with a positive definite advection scheme (Smolarkiewicz 
and Grabowski 1990). Results from the positive definite advection scheme are in better 
agreement with observations for tropical cloud systems (Johnson et al. 2002). 

In addition to the ARM surface fluxes, for land surface sensitivity experiments, the 
surface fluxes were extracted from the Land Information System (LIS, Kumar et al. 2004, 
http://lis.gsfc.nasa.gov/). LIS is a high performance land surface modeling and data 
assimilation system. It contains numerous land surface models (LSMs) that can be driven 
by a variety of atmospheric forcing from point to gridded data. For this study the NOAH 
LSM was employed. This LSM simulates soil moisture (both liquid and frozen), soil 
temperature, skin temperature, snowpack depth, snowpack water equivalent (and hence 
snowpack density), canopy water content, and the energy flux and water flux terms of the 
surface energy and surface water balances. The LSM land-surface parameters were 
initialized with University of Maryland 1 km datasets for vegetation and land-sea masks 
(Hansen et. al. 2005). Climatological datasets were ingested in order to initialize other 
vegetation parameters such as albedo and green vegetation fraction. Soils types were set 
using the State Soil Geographic Database for State [Soil Survey Staff, Natural Resources 
Conservation Service, United States Department of Agriculture], which has a 1 km 
horizontal resolution. Initial soil water and temperature profiles were also assigned 
according to climatology. 
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The LSM was integrated for 15 years up through the study period. For the period 1985 
through 1996, NCEP reanalysis data (the NCEP reanalysis data was obtained from the 
NOAA-CIRES ESRL/PSD Climate Diagnostics branch, Boulder, Colorado, USA, web 
site at http://Www.cdc.noaa.gov/) was used for the atmospheric forcing. After this period, 
l/8th degree atmospheric forcing was provided by the North American Land Data 
Assimilation System (NLDAS, Cosgrove et. al. 2003), which incorporates high 
resolution GOES radiation and Stage IV precipitation fields into the NCEP ETA Data 
Assimilation System (EDAS). Modeled fluxes and temperature fields were then 
evaluated against a variety of surface station data and found to be in excellent agreement 
with observations. The modeled latent and sensible heat fluxes were then extracted for 
use in the GCE model. 
b. Data 

Two cases are studied in the present paper. The first one is the ARM Spring 2000 IOP, 
which was also used in the ARM CPM WG Case 4 Study (Xie et al. 2005; Xu et al. 
2005). T his dataset starts at 1730 UTC 1 March and ends at 0830 UTC 22 March 2000. 
The second case covers the period from 2030 UTC 25 May to 830 UTC 14 June 2002 for 
the same ARM domain, which overlaps IHOP 2002. The two cases represent springtime 
and summertime midlatitude clouds, respectively. 

The ARM observational data used are classified into two parts: forcing and evaluation 
data. Large-scale forcing data (i.e., vertical motion and horizontal advective tendencies of 
temperature and moisture) are derived using data collected from the two ARM IOPs and 
the variational analysis approach described in Zhang and Lin (1997) and Zhang et al. 
(2001). The values represent the mean ARM Cloud and Radiation Testbed (CART) 
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domain rather than a single point (Zhang et al. 2001). The surface fluxes are obtained 
from site-wide averages of observed fluxes from the ARM Energy Balance Bowen Ratio 
(EBBR) stations. The fluxes are assumed to be horizontally uniform in the model. The 
LIS fluxes, which provide an alternate source for surface fluxes, are used for comparison 
and are discussed in section 4.b as well as in the preceding subsection. 

Evaluation data include observed temperature and humidity as well as others on 
clouds and precipitation. Temperature and humidity are observed every three hours 
during the IOPs. Cloud liquid water content and ice water content are obtained as ARM 
MICRO-BASE products (Miller et al. 2003). Vertical profiles of cloud fraction are 
derived from the hydrometeor frequencies from the Active Remotely Sensed Cloud 
Layers (ARSCL) data archive (Clothiaux et al. 2000). The uncertainties in those 
quantities were discussed by Xie et al. (2005). 
c. Numerical experiments 

A default numerical experiment is set up with a 1 km horizontal resolution, vertical 
resolution that ranges from 42.5 m at the bottom to 1 km at the top, and an integration 
time step of 6 s. The GCE model uses 128x128x41 gridpoints and is integrated for 20 
days except for specific tests. The model domain is located over the ARM Southern Great 
Plains (SGP) site with the center at 36. 6N and 96. 5W. 

The numerical experiments discussed in the paper are listed in Table 1. There are two 
control experiments COO and C02 for the 2000 and 2002 cases, respectively. Experiments 
COO, BOO and GOO are used to examine the sensitivity of clouds to domain structure. COO 
is a springtime cloud simulation. Its results are compared with those of other CRMs as 
part of a model inter-comparison. In the experiment, surface fluxes come from the ARM 
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observations (Fig. 1). Experiment BOO uses 256x256x41 gridpoints and a 2 km horizontal 
resolution for a bigger domain to test the sensitivity of clouds to domain size. In contrast 
to BOO, experiment GOO uses 256x256x41 gridpoints but a 1 km horizontal resolution to 
test the sensitivity of clouds to grid size. 

Two experiments test the sensitivity of clouds to cloud microphysics. Experiment 
MOOICE takes account of cloud ice sedimentation, and MOORH uses the microphysics 
parameterization scheme of Rutledge and Hobbs (1984). In addition, a 2D experiment 
D00 is used to assess the sensitivity to dimensionality, Unking the 3D GCE model to 
other 2D CRMs in the comparison. 

Control experiment C02, contrasting springtime COO case, is a summertime cloud 
simulation. As with COO, its surface fluxes come from ARM observations and are shown 
in Fig. 2. As compared by Figs. 1 and 2, the 2002 case possesses strong surface latent 
heat flux both in amplitude and daily average. 

In contrast to C02, experiment L02 uses surface fluxes from LIS, a land data 
assimilation system, that are then averaged horizontally for comparison. A snap-shot of 
the spatial distribution of the LIS surface fluxes at day 1 is shown in Fig. 3. The 
horizontally-averaged surface fluxes are also shown in Fig. 2 in comparison with the 
ARM surface fluxes. As shown in the figure, the latent heat fluxes from ARM and LIS 
are very close while the sensible heat fluxes are different, especially from day 4 to 9. 
Experiment L02 tests the sensitivity of clouds to surface fluxes from different sources. 

3. The 2000 case for springtime clouds 

a. Control experiment (COO) 
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The control numerical experiment (COO) falls within the spring of 2000; Part of the 
large-scale forcing data over the CART domain is displayed in Fig. 4, including large- 
scale horizontal temperature advection, horizontal advection of water vapor mixing ratio, 
and large-scale vertical velocity. The surface latent and sensible heat fluxes are shown in 
Fig. 1. These forcing data start at 1730 UTC 1 March and end on 21 March 2000. Various 
synoptic systems pass through the domain in the following sequence: a synoptic 
cyclogenesis event (1-4 March), a cold front leg (5-8), an upper-level trough (9-11), non- 
precipitating clouds (12-15), a cold front with frontognesis (15-19), and stationary fronts 
( 20 - 22 ). 

The 3D GCE model with Lin et al. (1983) microphysics is used to simulate the 
twenty-day period. Modeled surface precipitation is compared with observations in Fig. 5 
with thin and thick lines, respectively. In general, model surface precipitation agrees well 
with observations. Differences in intensity exist at days 6.5, 12.5 and 14.5. The 
accumulated surface precipitation amount is 6.3% smaller than was observed. Both 
modeled and observed surface precipitation rates have similar probability distribution 
functions (PDFs). 

Time-pressure cross sections of retrieved liquid and ice water contents are displayed 
in Fig. 6. Since retrievals have been well tested on thin non-precipitating clouds but not 
on thick precipitating clouds (Dong and Mace 2003), contour lines in the figure indicate 
the relative magnitude of the water contents. However, the contour line dividing zero 
from non-zero water contents shows clearly the extent of the clouds. Such distributions of 
water contents are used to evaluate the modeling results. 


10 



The domain-averaged water contents in the model are displayed in Fig. 7, where the 
liquid and ice water contents are defined as the total mixing ratios of liquid and ice water 
species, respectively. As shown in Figs. 6 and 7, the distribution of liquid water content 
in the model is similar to that retrieved. However, the distribution of model ice water 
content is quite different from the retrieved. Model ice water, in contrast to the retrieved, 
persists above 265 hpa as a residue of modeled clouds. 

Time-pressure cross sections of observed relative humidity and cloud fraction are 
shown in Fig. 8 and cross sections of the same variables from the model in Fig. 9. A grid 
box in the model is defined as “cloudy” when the radar reflectivity dBZ e a;-35, where the 
radar reflectivity is estimated from the mixing ratios of all water species using the 
algorithm of Luo et al. (2003). The distributions, as shown in Figs. 8 and 9, are similar 
for the main precipitation events such as those at day 1, 6.5, 9.5, and 15. However, the 
distribution of observed cloud fraction has many more fine structures than does the 
model. Modeled relative humidity, as shown in Figs. 8 and 9, is larger than was observed 
in the upper troposphere. The modeled relative humidity near 265 hpa increases gradually 
with time. It is around 50-60%, implying that air there is saturated with respect to ice. 
The spuriously large relative humidity above 265 hpa is associated with too much ice 
there. 

Cloud fraction, temperature, and liquid and ice water contents from the model and 
observations are averaged over twenty days to show their mean profiles as a function of 
pressure (Fig. 10). Average liquid and ice water contents in the model are smaller and 
larger than those retrieved, respectively, although the difference between observations 
and the model is partly due to the way in which the values were obtained. 
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Cloud amount in the model, as shown in Fig. 10, is larger than observed. The modeled 


cloud amount is -20% more than was observed from 800 to 250 hpa. The model air 
temperatures at 100 and 265 hpa are 7.8K lower and 3.5K higher than observations, 
respectively. The simulated near-surface air temperature is 5K higher than the 
observations. The temperature differences between the model and observations can be 
explained by the over prediction of cloud residue in the upper troposphere. Since cloud 
ice in the upper troposphere emits longwave radiation into space and absorbs upward 
longwave radiation from the air and land surface below, too much cloud ice near 165 hpa 
decreases the temperature at -100 hpa and increases the temperature at -265 hpa. 
Meanwhile, too much cloud ice in the upper troposphere also increases the downward 
longwave radiation, which in turn contributes partly to the over prediction of air 
temperature near the surface. The magnitude of the temperature error is close to that in 
the other models for continental clouds (e.g., Xie et al. 2005, Xu et al. 2005) and larger 
than that for marine clouds (Johnson et al. 2002). 

Figure 1 1 displays the twenty-day average profiles of water vapor mixing ratio and 
relative humidity against pressure. The difference in the mixing ratios between the model 
and observations decreases from 1.1 g/kg near the surface to zero at the tropopause. In 
contrast, the difference in relative humidity generally increases with height, from nearly 
zero at the surface to 30% at 265 hpa. The spuriously high relative humidity in the upper 
troposphere is associated with the over prediction of cloud ice in the upper troposphere, 
since cloud ice increases (or decreases) through deposition (or sublimation) when air is 
saturated (or unsaturated) with respect to ice. 
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b. Sensitivity experiments 

Experiment GOO contains 256x256x41 gridpoints to test the sensitivity of simulations 
to domain size. This experiment shows (figure omitted) that the modeled surface 
precipitation agrees well with observations just as in COO. Slight surface precipitation 
fluctuations such as at day 1.5 are weakened. Liquid and ice water contents, relative 
humidity, and cloudiness resemble those shown in Figs. 7 and 9 for the control 
experiment COO. 

Another experiment, BOO, was done using 256x256x41 gridpoints and a horizontal 
grid size of 2 km. The modeled surface precipitation in this experiment agrees well with 
observations, too. Slight surface precipitation fluctuations such as at day 1.5 are almost 
gone. The relative humidity from 265 to 165 hpa is 10% smaller than in COO (or GOO) 
only from day 1.5 to 5.5. Even though the relative humidity there is closer to 
observations, it is still higher than was observed. In addition, the relative humidity from 
165 to 100 hpa is close to that in COO, much larger than observed. Corresponding to the 
difference between BOO and COO in relative humidity, the cloud amount from 265 to 165 
hpa in BOO is larger than in COO, implying that vertical circulations occur near 265-165 
hpa and dehumidify the atmosphere there. 

In addition to the preceding experiments on domain structure, two numerical 
experiments were done to determine the sensitivity of cloud residue to cloud 
microphysics. In one of them, the microphysical scheme of Rutledge and Hobbs (1984) is 
used, where the fall speed of graupel is smaller than that of hail in the scheme of Lin et 
al. (1983). The scheme is more suitable for the simulation of tropical clouds. However, 
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the over prediction of cloud residue is still present in the upper troposphere, just as in 
COO. 

Cloud ice sedimentation can be an important process (e.g., Wu et al. 1999; Hong et ah 
2004). Experiment M00ICE accounts for the sedimentation of cloud ice (Starr and Cox 
1985). After introducing the sedimentation of cloud ice, cloud residue is improved. The 
difference between the model and observations in terms of relative humidity, in contrast 
to Fig. 9, increases slowly with time near 165 hpa. The distribution of cloud ice in the 
upper troposphere contains breaks in a time-pressure cross section. However, modeled 
relative humidity and cloud ice are still higher than observed, which is consistent with the 
results of Wu et al. (1999) regarding the sensitivity of cloud residue to cloud ice 
sedimentation. 

c. Dimensionality experiment (D00) 

Current MMFs mainly employ 2D CRMs for computational economy (e.g., 
Grabowski 2001; Khairoutdinov and Randall 2001; Chem et al. 2005). Also, Xie et al. 
(2005) and Xu et al. (2005) used four 2D CRMs to simulate two short periods from the 
2000 case. Thus, the 2D version of the GCE model is used to simulate the case for twenty 
days to address dimensionality sensitivity. 

The numerical experiment, referred to here as experiment D00, uses 512x41 
gridpoints. Its surface precipitation (Fig. 12) shows rapid fluctuations at day 1, 6.5, 14.5 
and 17 in contrast to observations and the 3D model. However, after being averaged over 
a long period, the surface precipitation rate and its accumulated amount are close to 
observations. The final accumulated precipitation amount after the 20-day integration is 
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5.6% smaller than that observed. This 5.6% difference is slightly smaller than the 
corresponding difference of 6.3% for the 3D model. 

Time-pressure cross sections of relative humidity and ice water content for D00 are 
displayed in Fig. 13. When compared with Figs. 8 and 9, the figure shows that the 
relative humidity in the upper troposphere changes irregularly, which is quite different 
from the observations or from the 3D model. In addition, the distribution of ice water in 
the model is different from that retrieved. It also differs from the 3D model by having 
periods without cloud ice in the upper troposphere. Another 2D experiment with 128x41 
gridpoints shows the sensitivity of cloud properties to domain size in 2D. The modeled 
surface precipitation is similar to that in Fig. 12 except for stronger surface precipitation 
fluctuations. The modeled clouds have distributions similar to those in Fig. 13 except for 
many fine structures in the middle troposphere. The fine structures associated with 
relative humidity and cloud ice in the middle troposphere are associated with the rapid 
fluctuation in modeled surface precipitation. The 2D simulations, in contrast to the 3D, 
show that fluctuations in surface precipitation are sensitive to domain size. 

Figure 14 displays time-average values of relative humidity, cloud fraction, and the 
horizontal variance of vertical velocity and temperature against pressure for experiments 
COO and D00 as well as the observations. The horizontal variance of vertical velocity and 
temperature in the upper troposphere is larger in the 2D experiment. This shows that the 
local vertical circulations in the upper troposphere are stronger in the 2D experiment than 
in the 3D model. Furthermore, the relative humidity in the 2D experiment is closer to 
observations on average than is the 3D. Cloud fractions in both the 2D and 3D 
experiments are larger than was observed. 
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The differences between the 2D and 3D simulations are attributed to the energy 
transfer between scales (Moeng et al. 1996) and buoyancy damping, where the latter is 
important near the tropopause because the energy transfer is associated with nonlinear 
momentum terms (e.g., Lesieur 1990) and the nonlinear terms due to convective cells are 
not important near the tropopause. Consider a temperature perturbation in a 3D and a 2D 
inviscid dry model. The perturbation in the 3D dry model damps down with time so that 
the atmosphere approaches buoyancy equilibrium while the perturbation in the 2D dry 
model does not (Mapes 1993; Nilsson and Emanuel 1999; Sobel and Bretherton 2000; 
Raymond and Zeng 2000). Therefore, when convective clouds initiate strong vertical 
oscillations, these oscillations persist in the 2D model but not in the 3D model. This 
results in a larger temperature variance in the upper troposphere in the 2D simulation 1 . 
This fundamental difference between the 2D and 3D models regarding the oscillations 
can explain the rapid fluctuation of surface precipitation and little water in the upper 
‘ troposphere in the 2D simulations. 

Xie et al. (2005) and Xu et al. (2005) used four 2D CRMs to simulate two short 
periods from the same case. Comparisons show that the 2D GCE model simulates surface 
precipitation intensity better than the other 2D models, although modeled precipitation in 

all models is good. Just like the other models, the 2D GCE model over-predicted cloud 

') 

ice and relative humidity in the upper troposphere. Its vertical profile of relative 
humidity, compared to observations, is very similar to that from the Advanced Regional 
Prediction System at Langley Research Center (ARPS/LaRC, Xu et al. 2005). The 

1 The horizontal variance of temperature and vertical velocity in the upper troposphere 
may be influenced by numerical techniques such as the sponge layer. 
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temperature error versus height in the 2D GCE model is similar to the other models, 
decreasing from 2.5 K at 965 hpa to -4 K at 265 hpa. The ice content in the 2D GCE 
model is closer to observations than in the other models; however, the cloud fraction in 
the 2D GCE is still larger than observed though smaller than the other models. The 
difference in cloud fraction between the 2D GCE and other CRMs is mainly attributed to 
the different definitions of cloudy area in the models. 

4. The 2002 case for summertime clouds 

a. Control experiment (C02) 

In contrast to the springtime case in the preceding section, a summertime case is 
studied in this section. Figure 15 displays some of the large-scale forcing data for the 
2002 case, namely, large-scale horizontal temperature advection, horizontal advection of 
water vapor mixing ratio, and large-scale vertical velocity. Figure 2 displays the ARM 
surface latent and sensible heat fluxes versus time. All data start at 2030 UTC 25 May 
2002 and last for 20 days. Compared with Fig. 4, Fig. 15 shows that the 2002 case has 
weaker large-scale forcing (e.g., large-scale vertical velocity in the planetary boundary 
layer). Figs. 1 and 2 show that the 2002 case has strong surface latent heat flux both in 
amplitude and daily average. 

The control experiment for the 2002 case (C02) is a twenty-day simulation with the 
same parameters as COO. Figure 16 displays the observed and modeled surface 
precipitation versus time. The strong precipitation events are fairly well captured. Just 
like in the other CRMs (Xu et al. 2002), the precipitation events at day 1.5 and 17 are 
delayed. To test the reason for the delay in precipitation, a numerical experiment was 
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made that started at day 1 with an artificial increase in water vapor in the planetary 
boundary layer. The results (figure omitted) show that the three precipitation events from 
day 1-3 are modeled well and the accumulated precipitation amount also agrees well with 
observations. This infers that the delay in precipitation on day 1.5 in C02 can be 
attributed to the lack of proper triggers for convective clouds. 

Two spurious precipitation events, as shown in Fig. 16, appear at day 13.5 and 14.2 in 
contrast to observations. The precipitation event at day 14.2 is initiated by large CAPE 
(convective available potential energy, see Fig. 17) and further intensified by the upward 
large-scale motion in the middle troposphere (see Fig. 15). To test the influence of 
previous accumulative errors on the spurious precipitation events, a new simulation was 
done that started at day 13. In this experiment (figure omitted), the precipitation event at 
day 13.5 disappears, and the precipitation event at day 14.2 is significantly weakened but 
still there. In summary, experiment C02 and other two experiments show that convective 
initiation, water spin up, and large-scale forcing influence the simulation of summertime 
precipitation. 

Figure 16 also shows that the accumulated precipitation amount in the model is 
smaller than the observations by 10.2%. Also, the PDF for the modeled precipitation rate 
is similar to observations at small rainfall rates but different at high rainfall rates. In 
general, the model has a reasonable accumulated precipitation and rainrate PDF 
compared to observations. 

Figure 17 displays time-pressure cross sections of observed relative humidity and 
cloud fraction, showing the diurnal variation of relative humidity and cloud amount in the 
lower troposphere. Figure 18 displays the same variables from the model. Based on these 
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two figures, the model relative humidity is larger than observations, especially in the 
upper troposphere. A pronounced diurnal variation of low clouds exists in the model. 

Twenty-day average values of water vapor mixing ratio and relative humidity (Fig. 
19) shows that observations and the model are similar in terms of water vapor mixing 
ratio but different for relative humidity. The modeled relative humidity is -30% higher 
than observed in the upper troposphere. Compared with Fig. 11, it shows that the 
differences between the model and observations for humidity are similar in both the 2000 
and 2002 cases although the surface water vapor mixing ratio and latent heat flux are 
much larger for the 2002 case. 

b. Land surface flux experiment (L02) 

Surface fluxes are necessary in an MMF to help drive large-scale circulations over 
continents. Unlike the numerical experiments in the preceding sections, no observational 
fluxes are available for MMF. Thus, it is of interest to test the sensitivity of cloud 
properties to surface fluxes when the surface fluxes are provided from a land surface 
model. This section describes such a sensitivity experiment and its impact on the 
simulation of clouds and precipitation. 

Experiment L02 follows the control experiment C02 except that surface fluxes as well 
as land surface temperature come from LIS. The LIS data are obtained after a 15-year 
spin-up driven by observational data with a 1 km horizontal resolution (see section 2. a for 
details). Figure 2 displays the domain- average values of the surface fluxes. The LIS 
fluxes have a strong diurnal signature very similar to the ARM data. The twenty-day 
average values of latent and sensible heat fluxes from LIS are 117.1 and 44.3 W/m 2 , 
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while the corresponding ARM data are 1 17.4 and 34.9 W/m 2 , respectively. However, LIS 
sensible heat fluxes from day 4 to 9 are larger than ARM ones in both amplitude and 
daily average. 

In experiment L02, LIS fluxes are assumed to be horizontally uniform. Surface 
precipitation characteristics (Fig. 20) are very similar to those in C02, but the rainrate 
PDF between the model and observations is different at high rainfall rate. Time-pressure 
cross sections of relative humidity and cloud amount from L02 are shown in Fig. 21. This 
and others figures on water contents (figure omitted) indicate that cloud residue is still 
over-predicted in the upper troposphere. However, the diurnal variation of cloud amount 
in the lower troposphere from day 4 to 9 is superior to C02. This improvement is 
attributed to the large LIS sensible heat flux from day 4 to 9 (see Fig. 2 for details). When 
the sensible heat flux increases, the upward flux of heat in the planetary boundary layer 
increases correspondingly, leading to the increase of the upward flux of water vapor 
there. Thus the surface relative humidity decreases due to the increase of air surface 
temperature and the decreases of water content near the surface, which increases the 
lifting condensation level and in turn decreases cloud amount in the lower troposphere. 
This connection between surface sensible heat flux and the diurnal variation of clouds in 
the lower troposphere (as shown in Figs. 18 and 21) is consistent with the difference in 
relative humidity in the planetary boundary layer between the two experiments. 

The ARM and LIS fluxes come from different sources. The ARM surface fluxes used 
in the control experiment come from the EBBR stations, which use the Bowen ratio to 
partition the fluxes. There are a total of 14 EBBR stations. A grid of 0.5 degree x 0.5 
degree was first set up to cover the ARM SCM (single-column model) domain. Next, the 
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Bames scheme with a length scale of 80 km was used to fill all the boxes. These 0.5 
degree x 0.5 degree boxes within the SCM domain were then averaged to get the area- 
averaged surface fluxes. Based on this procedure for the ARM surface flux data, it is 
inferred that some uncertainty is introduced into the ARM-averaged fluxes due to the 
small representative scale of land variables. In contrast, the LIS fluxes were obtained 
from a land surface model that was driven with observational data of 1 km-resolution (see 
section 2. a for details). Although the LIS fluxes at a point are not as accurate as the ARM 
fluxes, their fine resolution may lead to better area-averaged fluxes than those from 
“sparse” observational stations. However, this topic remains open and further 
comparisons between the two kinds of surface flux data are needed. 

5. Discussions 

Thermodynamics dominate atmospheric circulations on sufficiently large scales 
(Neelin and Held 1987; Emanuel 1995, 1999; Raymond 1995, 2000; Raymond and Zeng 
2000; Zeng et al. 2005). Consequently, large-scale vertical circulations are sensitive to 
atmospheric radiative cooling rate and the surface fluxes from the underlying surface 
(e.g., Zeng et al. 2005). Due to the significant influence of clouds on atmospheric 
radiation, the present paper compares clouds and surface fluxes in long-term simulations 
with observational data for the evaluation of CRMs in MMFs. 

Two twenty-day, continental midlatitude periods are simulated with the 3D GCE 
model and surface fluxes obtained from either observations or LIS. Observational large- 
scale forcing data from the ARM project are used to drive the model. Modeling results 
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are compared with the observed cloud data obtained from the ARM project, showing the 
benefits and bias of long-term CRM simulations. 

a. Cloud residue 

Modeled cloud residue (i.e., cloud water, ice, and water vapor) in long-term 
integrations is compared with observations from middle latitudes. Modeled cloud ice has 
a much higher frequency of occurrence than observations in the upper troposphere. 
Correspondingly, modeled relative humidity is much larger than observations there. A 
similar bias in cloud residue was evident in a 7-day simulation of GATE cloud systems 
(Grabowski et al. 1998). In addition, cloud fraction in the model is larger than 
observations. Modeled water and ice are smaller and larger than those retrieved, 
respectively, which is consistent with Xu et al. (2005) on the basis of 2D CRM 
simulations. 

Clouds, as shown in the sensitivity experiments, are influenced by many factors such 
as cloud microphysics, surface fluxes, dimensionality, and domain size. The over 
prediction of cloud residue may have multiple origins such as the periodic boundary 
condition as proposed by Grabowski et al. (1998) and Wu et al. (1998). 

b. Dimensionality 

Two-dimensional numerical experiments are performed to test the sensitivity of 
clouds and precipitation to dimensionality. As compared with observations, the 3D model 
produces better simulations than the 2D one. The 2D model has rapid fluctuations in 
surface precipitation. These fluctuations decrease when the number of domain gridpoints 
is increased. Grabowski et al. (1998) simulated oceanic clouds in the Tropics with both a 
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2D and 3D CRM and found the similar rapid fluctuations of precipitation in their 2D 
simulations in contrast to their 3D simulation. 

Relative humidity in the upper troposphere in the 2D model is higher than 
observations although the relative humidity there is closer to observations than in the 3D 
model. Modeled cloud ice in the upper troposphere covers a wider area than observations. 
The over prediction of cloud residue in the upper troposphere was also found in 2D 
simulations from four CRMs that were part of the ARM model inter-comparison (Xie et 
al. 2005; Xu et al. 2005). A similar over prediction of cloud residue occured in a 39-day 
2D simulation of TOGA COARE cloud systems (Wu et al. 1998). The vertical profile of 
relative humidity from the 2D GCE in comparison with observations is very similar to 
that from ARPS/LaRC (Xu et al. 2005). 

Clouds and precipitation are sensitive to dimensionality because of the difference in 
buoyancy damping and energy transfer between scales between 2D and 3D models. 
Temperature perturbations in the 3D model damp down due to gravity waves but not in 
the 2D model . Thus, oscillations in the 2D model, once generated by clouds, persist 
longer than in the 3D model, leading to cloud ice deposition in the upper troposphere and 
in turn dehumidification there after ice falls out. 

This evaluation of long-term CRM simulations provides a perspective on current 
MMFs. The over prediction of cloud residue (especially cloud ice) in the upper 

2 Readers interested in the difference in buoyancy damping between 2D and 3D models 
can start with a linear model for an inviscid and irrotational dry atmosphere. Once a 
temperature perturbation is imposed in the atmosphere, as shown by analytical solutions, 
it will have a different tendency in amplitude due to dimensionality. 
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troposphere, for example, partly explain “Red Spots” (or excessive precipitation) in the 
Tropics and over-active Madden- Julian oscillations in the current MMFs (e.g., 
Khairoutdinov and Randall 2001; Khairoutdinov et al. 2005; Chem et al. 2005), since 
clouds modulate radiation (e.g., Albrecht and Cox 1975; Grabowski et al. 1999) which in 
turn change large-scale circulations (e.g., Zeng et al. 2005) and consequently 
precipitation (e.g., Raymond and Zeng 2005) in a positive feedback (e.g., Raymond 2000; 
Raymond and Zeng 2000). 
c. Comparison with previous studies 

Several long-term CRM modeling studies have been performed recently, using the 
models with large-scale forcing. These previous studies can be compared to the present 
one to evaluate the sensitivity of cloud residue to model structure. 

Wu et al. (1998), using the Clark-Hall cloud model (Clark et al. 1996), performed 2D 
simulations of TOGA-COARE cloud systems for 39 days. They found that ice and liquid 
water contents were too large. Grabowski et al. (1998) noticed similar results in a 7-day 
simulation of GATE cloud systems from their 2D and 3D simulations. Wu et al. (1999) 
reported the sensitivity of cloud ice in the upper troposphere to cloud microphysics (or 
the sedimentation of cloud ice). They showed that the over prediction of relative humidity 
is improved with the increase of cloud ice sedimentation, but is still there. Similar results 
were obtained from other CRMs (Xu et al. 2002). The present results agree with the 
previous ones on the over prediction of cloud residue in the upper troposphere. 

Grabowski et al. (1998) and Khairoutdinov and Randall (2003) did long-term 
simulations, in contrast to observations, for dimensionality sensitivity. Grabowski et al. 
(1998) pointed out that 2D simulations produced a much higher temporal variability of 
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domain-averaged quantities in comparison to 3D simulations. The present simulations 
also support this modeling phenomenon. Moeng et al. (1996) showed a significant 
difference in horizontal variance of variables between the 2D and 3D simulations in the 
convective boundary layer. Khairoutdinov and Randall (2003) found that significant 
differences between 2D and 3D simulations had in temperature and vertical velocity 
variance, especially near the tropopause. They further inferred that a temperature 
perturbation in 2D does not decay with distance from its source as fast as in 3D. In the 
present study, the difference in vertical velocity and temperature variance between 2D 
and 3D simulations is also evident. These differences in variance, especially in the upper 
troposphere, can be explained by the difference in model buoyancy damping. 

It is not surprising that strong, local vertical circulations can influence cloud residue in 
the upper troposphere. Suppose that the air in the upper troposphere is saturated with 
respect to ice and that the ice falls out immediately as precipitation once deposition 
occurs. When strong vertical oscillations (or gravity waves) appear in the upper 
troposphere, the air is dehumidified, which shows how local vertical circulations can 
influence cloud residue. With this idealized case for perspective, it is not difficult to 
understand how differences in the present 2D and 3D simulations can impact cloud 
residue in the upper troposphere. 

6. Summary 

Two twenty-day, continental midlatitude cases are simulated with a CRM and 
compared with ARM data to evaluate the benefit and bias of CRMs in MMFs. Modeled 
surface precipitation is compared with observations first and then cloud residue is 
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compared, because surface precipitation and cloud residue are associated with cloud 
growth and decay, respectively, to some extent. 

All numerical experiments in the present paper output reasonable surface 
precipitation, as have other CRMs (e.g., Grabowski et al. 1998; Khairoutdinov and 
Randall 2003; Xie et al. 2002, 2005; Xu et al. 2002, 2005). Well-simulated surface 
precipitation implies that CRMs properly represent the upward flux of water vapor 
associated with cloud growth. However, all of the numerical experiments over-predict ice 
and relative humidity in the upper troposphere, which alters the vertical profile of 
temperature through atmospheric radiation. Such over prediction appears in other CRMs 
too (e.g., Xie et al. 2005; Xu et al. 2005) and needs to be addressed as a high priority for 
long-term CRM simulations. 

Special attention is given to the sensitivity of cloud properties to dimensionality in the 
present paper. Numerical results show that surface precipitation fluctuates more rapidly 
in the 2D CRM than in the 3D CRM. The 2D CRM also has less cloud residue in the 
upper troposphere. This sensitivity of simulated cloud properties can be explained by the 
difference in buoyancy damping between the 2D and the 3D models. 

Surface fluxes from LIS, a land data assimilation system, are compared with the ARM 
data. The LIS latent heat flux is close to observations, but its sensible heat flux is larger 
than observations for a period. When LIS surface flux data replace ARM data in the 
CRM simulations, similar results are obtained except that LIS brings about a better 
simulation of diurnal cloud variation in the lower troposphere. This work suggests that 
ARM and LIS surface flux data should be compared further with more cases in the future. 
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Caption 


Figure 1 Time series of surface fluxes for the 2000 case. The data start at 1730 UTC 1 
March 2000. Solid and dashed lines denote latent and sensible heat fluxes, respectively. 
Thick lines display daily-averaged values. 

Figure 2 Time series of surface fluxes for the 2002 case. All data start at 2030 UTC 25 
May 2002. Solid and dashed thin lines represent the surface fluxes from the ARM 
observations and LIS land data assimilation system, respectively. Thick lines represent 
corresponding daily-averaged values. 

Figure 3 Horizontal distribution of the surface fluxes at day 1 obtained from LIS. 

Figure 4 Time-pressure cross sections of large-scale horizontal temperature advection 
(upper), horizontal advection of water vapor mixing ratio (middle), and large-scale 
vertical velocity (lower pane). Data start at 1730 UTC 1 March 2000. Shaded areas 
indicate positive values; dashed and solid lines represent negative and positive valued 
contour levels, respectively. 

Figure 5 Surface precipitation for the spring 2000 period. Green and black solid lines 
represent observations and control experiment COO, respectively. Surface precipitation 
rate and accumulated rainfall are shown in the upper and middle panes, respectively. The 
PDF of surface precipitation is shown in the lower pane. 

Figure 6 Time-pressure cross sections of retrieved liquid and ice water contents starting 
from 1730 UTC 1 March 2000. 

Figure 7 Time-pressure cross sections of liquid and ice water contents obtained from the 
control experiment COO for the spring 2000 period. 
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Figure 8 Time-pressure cross sections of observed relative humidity and cloud fraction 
starting from 1730 UTC 1 March 2000. In the upper part, shaded areas indicate a relative 
humidity of more than 50%; dashed and solid lines represent the contour levels with 
relative humidity smaller and larger than 50%, respectively. 

Figure 9 Time-pressure cross sections of relative humidity and cloud fraction from the 
control experiment COO for the spring 2000 period. In the upper part, shaded areas 
indicate a relative humidity of more than 50%; dashed and solid lines represent the 
contour levels with relative humidity smaller and larger than 50%, respectively. 

Figure 10 Twenty-day mean profiles of liquid and ice water content, cloud fraction, and 
air temperature difference between the model and observations against pressure in the 
2000 case. Thick and thin lines represent variables from observations and experiment 
COO, respectively. 

Figure 11 Same as Fig. 10 except for relative humidity and the mixing ratio of water 
vapor. 

Figure 12 Same as Fig. 5 except for the experiment D00 (a 2D numerical experiment for 
the 2000 case). 

Figure 13 Time-pressure cross-sections of relative humidity and ice water content in 
experiment D00 (a 2D numerical experiment for the 2000 case). 

Figure 14 Twenty-day mean profiles of relative humidity, cloud fraction, vertical 
velocity variance, and temperature variance against pressure in the 2000 case. Thick lines 
represent observed variables; thin lines variables from experiment COO; and thin dashed 
lines variables from experiment D00. 


38 



Figure 15 Same as Fig. 4 except for the 2002 case. Data start at 2030 UTC 25 May 2002. 


Figure 16 Same as Fig. 5 except for the 2002 case. The numerical experiment is C02. 
Figure 17 Same as Fig. 8 except for the 2002 case. 

Figure 18 Same as Fig. 9 except for the 2002 case. Data are from the control experiment 
C02. 

Figure 19 Same as Fig. 1 1 except for the 2002 case. The numerical experiment is C02. 
Figure 20 Same as Fig. 16 except for the 2002 case and experiment L02. 

Figure 21 Same as Fig. 18 except for the 2002 case and experiment L02. 
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Table 1 List of Numerical Experiments 


Experiment 

Period 

Dimen 

Grid 

Gridpoint 

Surface 

Comments 

Name 


-sion 

Size 


Flux 


COO 

3/1-3/20,2000 

3 

1 km 

128x128x41 

ARM 

control run for case 
2000 

BOO 

3/1-3/20, 2000 

3 

2 km 

256x256x41 

ARM 

big domain 

GOO 

3/1-3/20, 2000 

3 

1 km 

256x256x41 

ARM 

spatial grid size 

MOOICE 

3/1-3/20, 2000 

3 

1 km 

128x128x41 

ARM 

run with 
sedimentation of 
cloud ice 

MOORH 

3/1-3/20, 2000 

3 

1 km 

128x128x41 

. 

ARM 

run with the scheme 
of Rutledge and 
Hobbs (1984) 

DOO 

3/1-3/20, 2000 

2 

1 km 

512x41 

ARM 

dimensionality 

sensitivity 

C02 

5/25-6/13, 2002 

3 

1 km 

128x128x41 

ARM 

control run for case 
2002 

L02 

5/25-6/13, 2002 

3 

1 km 

128x128x41 

LIS 

land surface flux 
sensitivity 


40 




























































Mar 1 - 20, 2000 



Time (day) 


Figure 1 Time series of surface fluxes for the 2000 case. The data start at 1730 UTC 1 
March 2000. Solid and dashed lines denote latent and sensible heat fluxes, respectively. 
Thick lines display daily- averaged values. 
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Figure 2 Time series of surface fluxes for the 2002 case. All data start at 2030 UTC 25 
May 2002. Solid and dashed thin lines represent the surface fluxes from the ARM 
observations and LIS land data assimilation system, respectively. Thick lines represent 
corresponding daily-averaged values. 
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Figure 3 Horizontal distribution of the surface fluxes at day 1 obtained from LIS. 
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Figure 4 Time-pressure cross sections of large-scale horizontal temperature advection 
(upper), horizontal advection of water vapor mixing ratio (middle), and large-scale 
vertical velocity (lower pane). Data start at 1730 UTC 1 March 2000. Shaded areas 
indicate positive values; dashed and solid lines represent negative and positive valued 
contour levels, respectively. 
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Figure 5 Surface precipitation for the spring 2000 period. Green and black solid lines 
represent observations and control experiment COO, respectively. Surface precipitation 
rate and accumulated rainfall are shown in the upper and middle panes, respectively. The 
PDF of surface precipitation is shown in the lower pane. 
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Figure 6 Time-pressure cross sections of retrieved liquid and ice water contents starting 
from 1730 UTC 1 March 2000. 
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Figure 8 Time-pressure cross sections of observed relative humidity and cloud fraction 
starting from 1730 UTC 1 March 2000. In the upper part, shaded areas indicate a relative 
humidity of more than 50%; dashed and solid lines represent the contour levels with 
relative humidity smaller and larger than 50%, respectively. 
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Figure 9 Time-pressure cross sections of relative humidity and cloud fraction from the 
control experiment COO for the spring 2000 period. In the upper part, shaded areas 
indicate a relative humidity of more than 50%; dashed and solid lines represent the 
contour levels with relative humidity smaller and larger than 50%, respectively. 
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Figure 10 Twenty-day mean profiles of liquid and ice water content, cloud fraction, and 
air temperature difference between the model and observations against pressure in the 
2000 case. Thick and thin lines represent variables from observations and experiment 
COO, respectively. 
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Figure 11 Same as Fig. 10 except for relative humidity and the mixing ratio of water 
vapor. 
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Figure 13 Time-pressure cross-sections of relative humidity and ice water content in 
experiment D00 (a 2D numerical experiment for the 2000 case). 
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Figure 14 Twenty-day mean profiles of relative humidity, cloud fraction, vertical 
velocity variance, and temperature variance against pressure in the 2000 case. Thick lines 
represent observed variables; thin lines variables from experiment COO; and thin dashed 
lines variables from experiment D00. 
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Figure 15 Same as Fig. 4 except for the 2002 case. Data start at 2030 UTC 25 May 2002. 


55 







Pressure (hpa) Pressure (hpa) 



{Obs) 


Contour interval = 10 


May 25- Jun 13, 2002 


Time (day) 


Relative Humidity (%) 
65 
165 
265 
365 
465 
565 
665 
765 
865 
965 



Pressure (hpa) Pressure (hpa) 



(Exp C02) 


Relative Humidity (%) 


65 

165 

265 

365 

465 

565 

665 

765 

865 

965 


Contour Interval = 10 


May25-Jun 13,2002 


Time (day) 



0 10 20 


May25-Jun 13,2002 


Time (day) 

BO 90 


65 

165 

265 

365 

465 

565 

665 

765 

865 

965 


Cloud Amount (%) 


(Exp C02) 


Figure 18 Same as Fig. 9 except for the 2002 case. Data are from the control experiment 
C02. 
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Figure 21 Same as Fig. 1 8 except for the 2002 case and experiment L02. 
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